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ABSTRACT 

The spin-down of a neutron star, e.g. due to magneto-dipole losses, results in compres¬ 
sion of the stellar matter and induces nuclear reactions at phase transitions between 
different nuclear species in the crust. We show that this mechanism is effective in 
heating recycled pulsars, in which the previous accretion process has already been 
compressing the crust, so it is not in nuclear equilibrium. We calculate the corre¬ 
sponding emissivity and confront it with available observations, showing that it might 
account for the likely thermal ultraviolet emission of PSR J0437—4715. 

Key words: stars: neutron - stars: interiors - pulsars - nuclear reactions ultraviolet: 
stars 


1 INTRODUCTION 

The detection of likely thermal ultraviolet emis¬ 
sion from the millisecond pulsar PSR J0437—4715 
(Kargaltsev, Pavlov & Romani 2004; Durant et al. 2012) 
indicates that some reheating mechanisms operate in old 
pulsars. Gonzalez & Reisenegger (2010) reviewed various 
possibilities, concluding that the most effective ones are 
vortex creep (Alpar, Pines, Anderson & Shaham 1984) 
and rotochemical heating (Reisenegger 1995), both very 
sensitive to nucleon superfluidity parameters. Although 
the theoretical and observational uncertainties are difficult 
to quantify (Gonzalez-Caniulef & Reisenegger, in prepa¬ 
ration), calculations by Petrovich & Reisenegger (2010) 
and Gonzalez-Jimenez, Petrovich & Reisenegger (2015) 
suggest that rotochemical heating can explain the observed 
temperature of PSR J0437—4715 only if either the proton 
or neutron energy gap is sufficiently large throughout the 
whole stellar core, which appears to be unlikely, at least 
for massive neutron stars (NSs). Here we discuss a variant 
of rotochemical heating that operates in the crusts of 
pulsars that were previously recycled through accretion. 
In these, the accreted matter has been compressed slowly, 
at moderate temperatures, and therefore its nuclear trans¬ 
formations have not been able to reach its ground state of 
fully catalyzed matter (Haensel & Zdunik 1990b). Thus, a 
small additional compression, now due to the decreasing 
centrifugal force, can immediately induce further reactions 
at the interfaces between layers of different dominant 


nuclei 1 . This process provides an additional heating source 
that has been neglected in the literature and is considered 
here. We emphasize that the proposed heating mechanism 
is not relevant for fully catalyzed (equilibrium) NS crusts, 
which need to be compressed substantially in order to ini¬ 
tiate these transformations. This explains why Iida & Sato 
(1997) found a much lower heating rate than that obtained 
in the present paper (see also footnote 2). 

In Section 2 we discuss how the compression produced 
as the star spins down initiates nuclear reactions at isobaric 
and isopycnic phase transition surfaces in the crust, leading 
to energy release and heating. In Section 3, we present three 
ways of calculating the heating rate, at different levels of ap¬ 
proximation: a simple, Newtonian estimate, a slow-rotation 
approximation within General Relativity, and a general- 
relativistic, numerical calculation valid also at fast rotation 
rates. Section 4 compares the three approaches, presents the 
results and analyzes the role of the proposed new heating 
mechanism in the interpretation of observations of the pul¬ 
sar PSR J0437—4715. We conclude in Section 5. 


1 The very similar process called ‘deep crustal heating’ occurs 
during the accretion process, in which the matter falling onto the 
NS surface compresses the underlying layers of the crust; see, e.g., 
Haensel & Zdunik 1990a; Brown, Bildsten & Rutledge 1998. 
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2 CRUST COMPRESSION AND HEATING IN 

A SPINNING-DOWN NS 

As a NS spins down, the compression and shear components 
of its deformation tensor are of the same order of magni¬ 
tude, but the shear modulus of the crust is much smaller 
than the bulk modulus (e.g., Haensel, Potekhin & Yakovlev 
2007), allowing us to neglect shear stresses and treat the 
whole star as a fluid. In the rotating star, the condition of 
hydrostatic equilibrium requires isobaric (constant pressure 
P ) and isopycnic (constant density p) surfaces to coincide 
(see, e.g., Fernandez & Reisenegger 2005). Since the relation 
between pressure and density depends on chemical compo¬ 
sition, the latter must also be uniform on each of these sur¬ 
faces, but can change from surface to surface. Thus, as the 
star spins down, any such surface will deform (becoming 
more spherical) and the associated pressure and density will 
increase, but the same particles must remain associated with 
it (otherwise the chemical composition would become non- 
uniform on a given isobaric surface as an NS decelerates) 2 . 

In the crust, there are various phase transitions (which 
we label by an index i), defined by critical pressures Pi. 
When the pressure of a given isobaric surface reaches Pi, 
nuclei at this surface undergo exothermic reactions (electron 
captures, neutron emissions, or pycnonuclear reactions), re¬ 
leasing an energy qi per nucleon. If the total number of nu¬ 
cleons at P < Pi is A Ni, and the angular velocity f 1 of the 
star decreases, the number of nucleons crossing the phase 
transition i per unit time is —(dANi/dQ)Cl, where we use 
the convention that total derivatives correspond to stellar 
models with the same total baryon number A (but gener¬ 
ally different rotation rates) and dots denote time deriva¬ 
tives. The total energy released per unit time at all these 
phase transitions, redshifted to a distant observer, is 

tpoo dANi Clqt f ^ 

E (1) 

i (0 

The gravitational redshift correction is the inverse of the 
time component of the four-velocity u M of a fluid element 
(Miralles, van Riper & Lattimer 1993), 

u = (- g t t - - O. 2 g^) 1/2 (2) 

where gtt, gt 4 >, and g are components of the metric of 
the rotating star. It can be shown that u 4 is also uniform 
on isobaric surfaces (Fernandez & Reisenegger 2005), so its 
value at a given transition i is well-defined. 


3 EVALUATING THE HEATING RATE 
3.1 A simple Newtonian estimate 


Applying the Newtonian condition of hydrostatic equilib¬ 
rium to the (assumed) thin layer of matter above the tran¬ 
sition i in a non-rotating star, we can write 

4n PiR 2 4t rRR 4 


A Ni{Q = 0) 


m B g GMmti ’ 


(3) 


2 This contradicts a result of Iida & Sato (1997), who found 
that the stellar matter compresses at the equator and expands 
at the poles. This suggests that there is some mistake in their 
approach, most probably in their calculations of the Lagrangian 
displacements (see their equation 52). 


where R, M, and g are the stellar radius, mass, and surface 
gravity, G is the gravitational constant, and m-B is the mass 
of a free nucleon (taken to be the same for neutrons and 
protons). When the star rotates slowly (f1 11 k, where 

17k ~ (GM/ i ? 3 ) 1 / 2 is the Kepler frequency or break-up spin 
rate), there will be a correction: 


A Ni(Q) 


AJVj(O) 



(4) 


with a being a positive coefficient of the order of unity. Thus, 

dANi ■ „ I7f7 . Ar 8ttR 7 ■ . . 

dfT n ~ AiVi(0) * G 2 A7 2 m B nnPi ’ (5) 

from which we get the approximate expression 

8ivR 7 


E « —a 


G 2 M 2 m B 


Pi dil 


( 6 ) 


where we also suppressed the redshift factor because we ig¬ 
nored the effects of general relativity. Also, for simplicity, we 
assumed the same value of a for all phase transitions. Scaling 
with fiducial values Re = R/(10 6 cm), A/ 1.4 = M / (IAMq), 
H 3 = f 2 /( 10 3 s _1 ), and fl_i 4 = S7/(10 -14 s -2 ), equation ( 6 ) 
becomes 


E 


-7 x 10 28 a 


R< 


ML 


173 17-14 


E 


Pi 


Qi 


10 31 erg cm -3 1 MeV 


ergs 


(7) 


Notice the strong dependence of the heating rate (and thus 
the bolometric luminosity) on the stellar radius (oc R 7 ). This 
is reduced, however, when considering the effective temper¬ 
ature, T eff oc ( E/R 2 ) 1/4 oc R s/4 , or the Rayleigh-Jeans flux 
(relevant in the ultraviolet region), Frj oc R 2 T e g oc R 13 ^ 4 . 

It follows from equation (7) that the heating rate gener¬ 
ally depends only on 4 “microphysics-related” parameters, 
M, R, Pi, and qt (and this is verified by the results of nu¬ 
merical analysis in Section 4). The first two (M and R) 
are mainly determined by the core equation of state (EOS), 
while the second pair (Pi and qi) depends on the (rather 
uncertain) details of the accreted crust EOS (see Section 4 
for a more detailed discussion). 


3.2 General-relativistic slow-rotation 
appr oximat ion 

Still in the limit of slow rotation, I7/I7k < 1, we now aim 
at a more precise evaluation including the effects of General 
Relativity. When an NS decelerates (f7 decreases), its cen¬ 
tral density p c , given by p c (S7) = p c o + <5p c (I7), increases (p c 0 
is the central density at 17 = 0 ; 5p c is the small f7-dependent 
correction). The correction 5p c can be found from the condi¬ 
tion that the total number of baryons in the star, A(p c , 17), 
remains unchanged as the star spins down, 

A(p c0 + 5p c , 17) = A(p c0 , 0 ). ( 8 ) 

Expanding the left-hand side of ( 8 ) in a Taylor series, one 
arrives at the following quadratic equation for <5p c (0), 

/3 5p c + 7 8pi + a!7 2 « 0, (9) 

where f5 = dA/dp c , 7 = ( 1 / 2 )d 2 A/dp 2 , a = (l/2)d 2 A/dfl 2 
and all derivatives are taken at p = p c o and 17 = 0. To derive 
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Figure 1. Heating rate in the crust normalized to 10 29 ergs -1 as a function of NS mass for three BSk EOSs (left panel) and three PAL 
EOSs (right panel). Thick solid lines are calculated in the approximation of slow rotation, using formula (1) and (10). Thin dashed lines 
represent the approximation in equation (7) with the choice a = 1/2. Triangles show the heating rate, calculated for the BSk21 EOS 
with the full numerical procedure, described in Section 3.3. All curves are plotted for v = 100 Hz and S7 = —10 —14 s -2 . 


(9) we used the fact that dA/dSl oc 0 and hence vanishes at 
0 = 0. Knowledge of <5p c (0) (for more detailed analysis of 
equation 9 see below) allows us to write 

dANj _ dANj dp c dANj _ dANj d8p c dANj 

dSl dpc dfl + dfl dpc dSl dfl 

We evaluate the latter expression using the scheme of Hartle 
(1967) and Hartle & Thorne (1968), following the detailed 
description of Fernandez & Reisenegger (2005) and making 
use of our equation (9) and their equations (18)-(26), which 
we checked by an independent calculation. 

In the limit O —» 0, it is expected that dANi/dp c —> 
constant and dANi/dSl oc H, which we confirmed by nu¬ 
merical evaluations. On the other hand, dSpc/dSl has a 
less simple behaviour. As follows from equation (9), if /3 
is not small then the term 7 dp 2 can be neglected and 
Spc ~ —aSl 2 /P oc H 2 . As a result, E°° oc fill at S2 —» 0, 
as in equations ( 6 ) and (7). However, near the maximum 
(non-rotating) NS mass M max , which occurs at the same cen¬ 
tral density, p c , max, as the maximum baryon number, A max 
(e.g., Haensel, Potekhin & Yakovlev 2007), the situation is 
different. The coefficient /3 there is small (and vanishes at 
p c0 = p C!max ) so that the terms oc 5p c and 5p 2 in equation 
(9) can be of comparable magnitude. In the most extreme 
case (/3 = 0) one obtains Sp c = —{—aSi 2 / 7) 1 / 2 oc H so that 
E°° oc 12 at ^ 0 rather than oc fffl. 


3.3 General-relativistic, rapidly rotating NSs 

This problem was also studied numerically, without 
reference to a slow-rotation approximation, using the 
RNS code developed by S. Morsink and N. Stergioulas 

(http://www.gravity.phys.uwm.edu/rns/). This code 


constructs relativistic models of rapidly rotating NSs for a 
given tabulated equation of state. 

We calculate a set of NS models for the same total 
baryon number A and different spin rates Si. Then we look 
for the number of baryons located above the phase transition 
i making use of equation (20) of Fernandez & Reisenegger 
(2005), and approximate its dependence on SI by an analytic 
function. Differentiating this fit, we calculate dANi/dSl, to 
be used in equation ( 1 ). 

We follow the same strategy as in the slow-rotation limit 
(see Section 3.2), using equation (1) with u^y = (— gtt)^ 1 ^ 2 
evaluated for a non-rotating star. This approximation ap¬ 
pears to be rather accurate even for the most rapidly ro¬ 
tating NSs, not deviating from the accurate calculation by 
more than a couple per cent, as confirmed by Baubock et al. 
(2015). 


4 RESULTS 

To illustrate our results, we employ three EOSs of the BSk 
(Brussels-Skyrme) family (BSkl9, BSk20, and BSk21; see 
Potekhin, Fantina, Charnel, Pearson & Goriely 2013 for de¬ 
tails), and three stiff (compression modulus I\ = 240 MeV) 
EOSs from Prakash, Ainsworth & Lattimer (1988); here¬ 
after PAL EOSs (PAL I, PAL II, and PAL III) to describe 
matter in the NS core 3 . The unified BSk EOSs describe NS 
matter in the whole range of densities including the core and 
crust. However, Potekhin et al. (2013) do their calculations 

3 Softer EOSs are inconsistent with the precisely measured 
NS masses around 2Mq (Demorest et al. 2010; Antoniadis et al. 
2013). 
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only for catalyzed matter of NS crust, while here we are in¬ 
terested in NSs with an accreted crust that has not relaxed 
to the nuclear equilibrium state. Thus, although the actual 
functional dependence of P on p in the crust has little effect 
on our results (because the heating rate depends on M and 
R, see equation 7, which are mainly determined by the func¬ 
tion P(p) in the core and not in the crust), for consistency 
we employ an EOS of accreted crust from Haensel & Zdunik 
(1990b), which we match with the BSk or PAL EOSs in the 
core. On the other hand, the pressure at the phase tran¬ 
sitions (Pi) and the heat release per baryon qi in nuclear 
transformations at these transitions depend essentially on 
details of EOS in the accreted crust; they were taken from 
Table A.3 of Haensel & Zdunik (2008), thus assuming that 
X-ray ashes consist of pure 56 Fe. Note that, as is evident 
from the above discussion, the EOSs employed by us in the 
core and in the accreted crust have not been derived from the 
same microphysics input. Such calculations, which would al¬ 
low one to determine self-consistently the parameters M, R, 
Pi, and qt, are still unavailable in the literature (but one can 
make a rough guess about the possible uncertainty involved 
into the problem; see footnote 4 and a discussion below). 

Figure 1 allows to compare the results for all three ap¬ 
proaches described in the previous section. As it should be 
for a rather slowly rotating star (v = 100 Hz), the numerical 
results (triangles), calculated as discussed in Section 3.3, al¬ 
most coincide with the slow-rotation approximation (thick 
curves). Even the rough estimate from equation (7) (thin 
dashed curves) follows them quite well, only failing to re¬ 
produce the sharp upturn close to the maximum mass, due 
to the anomalous behaviour of the function dANi/dtt near 
the maximum mass discussed in Section 3.2, yielding just 
E°° oc Cl, in this limit, as opposed to E°° oc flCl at the lower 
masses (M < 0.95 M max ) to which the estimate from equa¬ 
tion (7) applies. Keeping these proportionalities in mind, it 
is easy to rescale this figure to any values of v -C Ok/( 27 t) 
and 12. 

As suggested by the estimate of equation (7) and con¬ 
firmed by the calculation results in figure 1, the heating rate 
E°° is indeed rather sensitive to the stellar mass and the 
EOS in the NS core (which together determine the stellar 
radius). It is also sensitive to the positions of phase transi¬ 
tions in the crust and energy release due to nuclear trans¬ 
formations at the phase transitions (Pi and qi, see Eq. 7). 
These quantities are rather model-dependent and, in partic¬ 
ular, depending on the nuclear symmetry energy, the heating 
rate E°° can vary by about a factor of two (Steiner 2012) 4 . 
Taking into account all these uncertainties, E°° can be esti¬ 
mated as 15°° (0.5 — 20) x lO^-Elot (where E IO t = IQCl is 
the rotation energy loss; I is the stellar moment of inertia). 

Figure 2 illustrates the accuracy of the slow-rotation ap¬ 
proximation as a function of rotation rate. As expected, for 
rapidly rotating stars (v ^ 500 Hz) the results of these two 


4 Actually, what Steiner (2012) estimated is the uncertainty of 
the deep crustal heating rate in accreting NSs, which is oc , qi. 
In our case, E°° is proportional to 54 , Piqi (see equation 7), so 
that, strictly speaking, the uncertainty can be somewhat differ¬ 
ent. However, Steiner (2012) does not provide any tables with his 
data for Pi and qi (as Haensel &; Zdunik 2008 did), so we cannot 
implement his results directly to check the real uncertainty. 


V, Hz 


0 200 400 600 800 





Figure 2. E°° normalized to 10 29 erg s —1 as a function of 12 
in units of Kepler frequency 12 k for BSk21 EOS and an NS with 
M = 1.75Mq (according to the RNS code 12 k = 10807.9 s —1 ). 
The thin curve is calculated in the slow-rotation approximation, 
using formulas (1) and (10), while the thick curve is calculated 
with the RNS code (see Section 3.3). The triangle shows E°° for 
v = 100 Hz, the frequency used in figure 1. To plot the figure we 
take 12 = —10” 14 s —2 . 

approaches differ substantially, with the exact solution being 
several-fold larger than the slow-rotation limit predicts. 

Surface temperatures of millisecond pulsars are gen¬ 
erally not known. However, Kargaltsev et al. (2004) and 
Durant et al. (2012) have measured likely thermal emission 
from PSR J0437—4715 in the far ultraviolet, implying that 
some reheating mechanisms must operate in this pulsar. The 
latter authors derive a redshifted effective temperature in 
the range 1.25 x 10 5 K < T^ < 3.5 x 10 5 K, under the 
assumption of a blackbody emitter with redshifted circum¬ 
ferential radius in the range 7.8 km < R 00 < 15 km. 

Figure 3 shows the predictions from the reheating mech¬ 
anism proposed in this paper. As one can see, rotation- 
induced heating in the crust is a powerful mechanism, be¬ 
ing close to explaining by itself (apart from other heating 
processes, such as vortex friction or rotochemical heating) 
the observed temperature of PSR J0437—4715. Note that 
the upper limit for R^, = 15 km assumed by Durant et al. 
(2012) when deriving the temperature from the ultravio¬ 
let flux is smaller than the values obtained for most of the 
currently considered EOSs, including those used in this pa¬ 
per. Since the ultraviolet flux (in the Rayleigh-Jeans limit 
of the thermal spectrum) is roughly oc R%aT((§, one could 
extrapolate the lower limit of their temperature range as 
T£ « 1.25x10 s (15km/R^fK, yielding «0.97xl0 s K 

at Roo = 17 km. Given various uncertainties in both the 
measurements and the theoretical models, as well as the 
crude assumption of a blackbody atmosphere, the predic¬ 
tions from our model might be considered to be in marginal 
agreement with the observation. The agreement might be 
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Figure 3. Effective TJg (in units of 10 s K) as a function 
of redshifted circumferential radius (or “radiation radius”) Roo 
for six EOSs (see figure 1) calculated in the slow-rotation ap¬ 
proximation with and H taken for PSR J0437—4715 (fl = 
1091s -1 , = 2.60 x 10 -15 s -2 ). Thin dashes show the mass 
range 1 M© < M < M max for each EOS, while thick lines 
correspond to 1.56M© < M < 1.96M© (the allowed mass 
range for PSR J0437—4715 at 68% confidence as measured by 
Verbiest et al. 2008). Arrows indicate the direction of increasing 
mass. 

improved by including some other reheating mechanisms as 
well. 


5 CONCLUSIONS 

We have shown that the compression of the accreted crust in 
the course of NS spin-down results in a significant energy re¬ 
lease at the phase transitions in the crust. The corresponding 
heating rate, E °°, can be estimated as (0.5 — 20) x 10 5 £ rot . 
In contrast to the other most effective reheating mecha¬ 
nisms proposed before — rotochemical heating (Reisenegger 
1995) and vortex creep (Alpar et al. 1984) — the reported 
reheating mechanism is not affected by the rather uncertain 
parameters of baryon superfluidity in the NS. As we have 
demonstrated, E°° mainly depends on the stellar mass and 
radius, as well as on the positions of the phase transitions 
in the crust and the energy release at these transitions; see 
equation (7). Noteworthy, it is insensitive to other details of 
the crust and core EOSs. 

The derived E°° is one to two orders of mag¬ 
nitude larger than the values obtained by Iida & Sato 
(1997) and is comparable to the rotochemical heating rate 
(Fernandez & Reisenegger 2005; Gonzalez-Jimenez et al. 
2015). For some EOSs, the predicted surface temper¬ 
ature is marginally consistent with the observation of 
PSR J0437—4715 (Durant et al. 2012). However, some other 
reheating mechanisms could also be important for this pul¬ 
sar. 
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